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We show here that the standard physical model used by Vriend et al. to analyse seismograph 
data, namely a non-dispersive bulk propagation, does not apply to the surface layer of sand dunes. 
According to several experimental, theoretical and field results, the only possible propagation of 
sound waves in a dry sand bed under gravity is through an infinite, yet discrete, number of dispersive 
surface modes. Besides, we present a series of evidences, most of which have already been published 
in the literature, that the frequency of booming avalanches is not controlled by any resonance as 
argued in this article. In particular, plotting the data provided by Vriend et al. as a table, it turns 
out that they do not present any correlation between the booming frequency and their estimate of 
the resonant frequency. 

PACS numbers: 



ACOUSTICS IN SAND DUNES 

It is a well-known fact that sand dunes present a lay- 
ered structure due to successive avalanches and capillary 
water retention. We focus here on the acoustic propa- 
gation within the metre scale surface layer composed of 
dry sand. According to Vriend et al. [2007], in [An- 
dreotti 2004; Bonneau et al. 2007], we "only consider 
low-speed surface waves of around 50 m/s" while "boom- 
ing results from the propagation of body waves, not sur- 
face waves" . This statement reflects a misunderstanding 
of acoustic waves in granular media. In fact, we have 
shown that, contrarily to ordinary elastic solids, bulk 
modes (the 'body waves') do not exist at all [Bonneau et 
al. 2007]. Instead, there is an infinite yet discrete num- 
ber of surface modes (not a single one) with a dispersive 
propagation (i.e. not a single wave speed independent of 
vibration frequencies). This has been directly evidenced 
in the field and in a lab experiment (Figure 1), and re- 
produced by an independent team [Jacob et al. 2007]. 

The explanation is simple. Sand is a divided medium 
that presents non-linearities of geometrical origin. As 
the grains do not have plane/plane contact, they behave 
like a spring whose stiffness depends on the normal force 
applied to put them into contact. As a consequence, 
the speed of sound c depends on confining pressure P - 
roughly as c oc Z 1 / 3 / 31 / 6 , where Z is the effective number 
of contacts per grain. A striking consequence is that the 
order of magnitude of the propagation speed c under a 
pressure of <~ 1 m of sand is lower than the sound veloc- 
ity in air (although density is 10 3 times larger). There- 
fore the effective compressibility is extraordinarily small! 
This property cannot be explained without involving ge- 
ometrical effects at the scale of a grain. This dependence 
on pressure has been proven several times experimentally 
and numerically for moderately large P (<~ 10 — 100 m 
of sand in Jia et al). Due to gravity-induced pressure 
(P oc pgz) , the surface layer of a dune does not constitute 
a homogeneous system as c increases with depth. Thus, 



no plane wave Fourier mode can exist in such a medium; 
only an infinite number of surface modes guided by the 
sound speed gradient may propagate. Consequently, a 
mode labelled n is localized within a depth n A below the 
surface (A is the wavelength) i.e. in a zone where the typ- 
ical pressure is P ~ pgnX. Thus, the surface propagation 
velocity increases with the mode number n in a similar 
way as c(P) (between n 1 / 4 and n 1 / 3 typically). So, even 
in the limit of an infinite depth (no finite-depth layering) , 
gravity produces a wave-guide effect, but no resonance. 

Considering now a layer of sand of finite depth H, a 
second wave-guide effect get superimposed to the first 
one and resonant modes may appear. As the system 
is still not homogeneous due to gravity, these resonant 
modes are not Fourier modes. By definition, they do not 
propagate - i.e. they have a vanishing group velocity - so 
that they correspond to cut-off frequencies of the system: 
no wave can propagate at a velocity smaller than the first 
resonant frequency Jr. The influence of the finite depth 
H is in fact limited to a very narrow range of frequencies. 
As soon as the depth H is larger than the wavelength 
(in practice, for a frequency / 25% above /#), it can 
be considered as infinite and the gravity effect prevails 
[Bonneau et al. 2007]. 

Then, it is easy to realize that if one strikes such a 
gradient-index medium with a sledge hammer (basically, 
the procedure used by Vriend et al.), a series of wave 
packets corresponding to the different surface modes will 
be propagated that are related to the gravity induced 
index gradient and not to the effect of the finite depth 
H. Moreover, this procedure only gives access to the 
group velocity at the mean frequency of excitation - yet 
not controlled nor specified by Vriend et al. - and not 
at all at the frequency / of spontaneous booming. Con- 
sequently, from such a procedure, just by reading the 
multi-modal structure of a seismograph obtained using 
only surface transducers, one cannot conclude to the ex- 
istence of multiple layers nor determine the resonant fre- 
quencies. Actually, if Vriend et al. were right, they could 
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easily provide a crucial test by performing measurements 
with transducers buried in the bulk, as those reported in 
Andreotti [2004] (see Auxiliary figures). The vibration 
in the soil during booming avalanches is already strongly 
reduced at 60 mm below the surface (Auxiliary Figure), 
which evidences directly that the song of dunes is not re- 
lated to a resonant mode at the meter scale but to surface 
propagation. 

As a conclusion, the quantity determined by Vriend 
et al. has probably nothing to do with the resonant 
frequency fn- Besides, one may wonder why they did 
not measured directly the resonant frequencies fn for the 
seek of comparison with the booming frequency /. Us- 
ing two different methods, we have ourselves performed 
such measurements using two different methods (Auxil- 
iary Material). For a depth H of dry sand of the order 
of 50 cm, we find typically 70 Hz; for 1 m, the resonance 
becomes hardly visible and is below 50 Hz. This is much 
smaller than the values found by Vriend et al. and than 
the booming frequency / (100 Hz in Morocco). So, as al- 
ready evidenced in our previous papers, the phenomenon 
is not driven by a resonance effect. As shown in Bonneau 
et al. [2007], the resonant frequency plays another role. 
As the surface waves do not propagate below fn, the later 
controls the threshold for the booming phenomenon. A 
large enough layer of dry sand is required for the surface 
waves to propagate and thus for booming to occur. The 
best situation is thus an infinite layer of dry sand and 
not a layered one. 

FIELD EVIDENCES AGAINST THE SELECTION 
OF THE FREQUENCY BY A RESONANCE 

Vriend et al. have missed another very important 
step of the argument presented in Andreotti [2004]. The 
booming frequency / is controlled by the shear rate 7 
(i.e. by the rate at which the grains collide with each 
other) in the shear zone separating the avalanche from 
the static part of the dune. Of course, 7 depends on the 
way the avalanche is forced. When a granular flow is 
driven by a pressure gradient, e.g. by pushing with the 
hand or the bottom, 7 can be varied and so does the emis- 
sion frequency /. Besides, the reproduction of the phe- 
nomenon at small scale [Haff 1986; Douady et al. 2006; 
Bonneau et al. 2007] is a direct proof that the acoustic 
emission is not related to the dune itself. Note in par- 
ticular that controlled lab experiments have allowed to 
produce sustained booming (not short squeaks), varying 
continuously from gravity to pressure gradient driving. 

Spontaneous avalanches on the slip face of dunes are 
driven by gravity in an homogeneous and steady way. 
In that case, only, 7 and thus / are expected to scale 
as \fgjd, with subdominant dependencies on the nature 
of the grains and the presence of cohesion. These con- 
ditions can be reproduced by a man-made slide only if 



the velocity of the body is constant. This is never the 
case in the movies provided by Vriend et al. to illus- 
trate their measurements. In the auxiliary animation, 
we show that a pulsed driving leads to a low frequency 
squeak similar to that measured by Vriend et al. (e.g. 
point at 70 Hz in their figure 3). It is pretty clear that 
such an inhomogeneous and unsteady driving do not re- 
semble the spontaneous avalanching process and thus, 
the avalanche has no reason to yield a constant booming 
frequency /. By contrast, during our 13 field trips in Mo- 
rocco from 2002 up to now, we have recorded around 20 
spontaneous avalanches and 100 man-triggered slides, on 
barkhan dunes ranging from 4 to 42 m in height. This was 
done in different places and by different weather condi- 
tions. We always measured the same frequency / within 
a tone. In particular, during our field trip of April 2007, 
we have performed a series of experiments for different 
flow thicknesses H, measured coherently by two indepen- 
dent techniques (see Auxiliary figures). Again, provided 
that avalanches were homogeneous and steady, the fre- 
quency / turned out to be constant. This demonstrates 
that the frequency is not either selected by a resonance 
effect over the depth of the avalanching flow. 



Figure 2 compares our predictions for the frequency 
/ of homogeneous booming avalanches (/ ~ 0.4-^/ g/d) 
and that proposed by the Caltech team (/ = fn). Once 
the data obtained by Vriend et al. in Dumont are plot- 
ted (fig. 2b) - and not anymore presented in a Table - 
no correlation can be observed between / and what they 
claim to be the resonant frequency fn (sec above). By 
contrast, with a choice of representation different from 
Vriend et al. (see tables and fig. 4 in Auxiliary Mate- 
rial), the relation between / and \J g/d appears beyond 
any doubt (fig. 2a and solid line in fig. 2b). This accu- 
mulation of evidences shows that the emission frequency 
/ is neither controlled by a resonance effect at the scale 
the avalanche as suggested by Douady et al. [2006] nor 
at the scale of a superficial layer as suggested by Vriend 
et al [2007]. 



The instability mechanism proposed in Andreotti 
[2004] is based on the interaction between the plastic 
deformations (shear band) and the elastic compression 
(acoustics) of the granular material. The collisions be- 
tween grains induce a transfer of momentum from trans- 
lation to vibration modes. In turn, the surface waves tend 
to synchronise the collisions. Importantly, this phase syn- 
chronisation linear instability do not depend on the fact 
that the vibrations modes are propagating or resonant. 
Although some 'mysteries' about singing dunes still re- 
main to be solved, our explanation remains the single one 
consistent with all the existing observations. 
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FIG. 1: Field experiment: propagation of a synthetic gaus- 
sian wave-packet at 350 Hz at the surface of a 50 cm deep 
fiat booming sand bed (proto-dune on a limestone plateau) 
i.e. without any layering effect, (a) Signal envelope mea- 
sured by correlation with the emitted signal at 5 cm (dotted 
line), 150 cm (solid line) from the source, (b) Comparison 
with the model proposed Bonneau et al. 2007. The three 
wave-packets received are direct evidences of the existence of 
multiple propagative surface modes. 



FIG. 2: (a) Booming frequency / as a function of grain size 
d. The plain symbols («)corresponds to homogeneous and 
steady avalanches (spontaneous or not) and the open ones 
(o) to the data of Vriend et al. (see details in the Auxiliary 
Tables), (b) Plot of the data-set presented by Vriend et al. as 
a Table (Large Dumont dune). Frequency / of man-induced 
avalanches as a function of the (biased) resonant frequency fn 
(see text). On both plots a) and b), the solid line corresponds 
to our prediction. The dotted line is the prediction of Vriend 
et al. [2007] 



AUXILIARY MATERIAL 

Auxiliary material for this comment gives the details 
necessary to understand our arguments and reproduce 
our measurements. There are three points of disagree- 
ment with Vriend et al. First, we show with a movie 
(Animation SI) that a man-made slide produces a boom- 
ing similar to spontaneous avalanches only if the body is 
kept rigid and the sliding velocity as constant as possi- 
ble. This movie is to be compared to those provided by 
Vriend et al. as Auxiliary Material, where the avalanches 
are systematically pulsed and inhomogeneous. Second, 
we show that the representation used by Vriend et al. 
to compare the frequency and the grain diameter (Fig. 
4) does not allow to see the obvious correlation of Fig. 
2a. We provide the table of the frequencies used in the 
figures of the comment (Fig. 2a, 2b and 4). A movie 
(Animation SI), a schematic of the measurement set-up 
(Fig. 5) , granulometric measurements (Fig. 3) and differ- 
ent frequency spectra (Fig. 6) are provided to document 
the data points obtained in the Atlantic Sahara. It may 
serve as a reference to start a data-base of steady homo- 
geneous booming avalanches. Third, we show that the 
method used by Vriend et al. to measure the first reso- 
nant frequency is not correct. We provide the measure- 
ment set-up for our controlled acoustic field experiments 
(Fig. 7a), including two methods to determine directly 
the resonant frequency (Fig. 8). Using these methods, 
we obtain a resonant frequency lower and not equal to 
the booming frequency. 
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Place 


f (Hz) 


d (fj, m) 


Kelso 


92.8 


200 


Sidi-Aghfenir 


103 


163 


Copiapo 


87 


210 


El Cerro Bramador 


75 


270 


Sand moutain 


63 


340 


Tarfaya 


90 


183 



TABLE I: Data points for which we do know that (i) 
avalanches where spontaneous or at least homogeneous and 
steady (ii) the grain diameter was determined from samples 
taken in the middle of the slip face, at the place where the 
avalanche was recorded. When possible, the frequency was 
averaged over several realisations. 



Place 


f (Hz) 


Big dune 
Eureka 
Kelso 
Dumont 


93.5 
79.3 
103.6 
86.4 



TABLE II: Average frequency / of the data points measured 
from Vriend et al. No description of the conditions of experi- 
ment has been provided by these authors. The six movies they 
provided show the same technique to create avalanches, which 
are pulsed and inhomogeneous. No description of the place 
where the grains were sampled has been provided by these au- 
thors. The data obtained in Dumont on December 2005 have 
not been taken into account for several reasons: no recording 
is provided, no movie is provided, the width of the frequency 
peak is extremely large (40% of the central frequency) which 
is the signature of a small squeak induced by an unsteady 
avalanche (see our Auxiliary movie). So it cannot be used 
in a plot aiming to show data obtained in homogeneous and 
steady conditions. 
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FIG. 3: Probability density function (PDF), weighted in 
mass, of the logarithm of the grain size d. Red dots: 4 m 
high barchan dune in Tarfaya that booms at the frequency 
/ = 90 Hz (see Auxiliary film). Black squares: 42m high 
mega-barkhan in Sidi Aghfinir that boomed at a mean fre- 
quency / = 103 Hz over 5 years (see Auxiliary film). The best 
fit by a log-normal distribution (a parabola in that represen- 
tation) gives the centre of the distribution: d — 183 ± 2 /im 
and d — 163 ± 2 /im respectively. Inset: Photograph of few 
grains from the small barchan dune in Tarfaya. 
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FIG. 4: How can the same d ata s how a clear scaling of the 
booming frequency / with \fgjd and, as stated by Vriend 
et al, no correlation? (a) Replot of figure 2a. The prediction 
corresponds to a straight line at 45° in the graph (equal range 
of axis). The symbols show the conditions of booming: ho- 
mogeneous and steady with a grain diameter d determined in 
the middle of the avalanche slip face (•) or pulsed and an un- 
known position at which sand was sampled (o). (b) Same but 
plotting d as a function of /. A large axis range is chosen for 
d (0 to 600 (im while the dune grain size ranges from 150 fim 
to 350 £tm), but the same narrow axis range for /. The solid 
line shows the prediction, (c) Same but adding error bars to 
indicate the width of the grain size distribution but letting 
no error bars on /, removing the last point on the left (Sand 
mountain), which significantly contributes to extend the dy- 
namics, and removing the Drediction. One obtains the Dlot 
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FIG. 5: (a) Schematic of the measurement set-up used to 
characterise booming avalanches. The pressure in the air is 
measured using a calibrated pressure transducer at 10 cm 
above the soil. A sound level apparatus using another mi- 
crophone is filmed by the CCD camera (see Animation SI). 
The vibration of the grains is measured with piezo-electric 
accelerometers that can be put at different depth. The three 
measurements of the booming amplitude are checked to be 
consistent with the loud-speaker model (Andreotti 2004). The 
flowing depth is measured with two independent techniques. 
First, we bury a piece of string (green on the schematic) of 
diameter 4 mm that can be curved with a very small force, 
attached at one hand to a large heavy object. The avalanch- 
ing grains entrain the string and align it along the steepest 
slope, at the surface of the static part of the dune. Second, we 
prepare in advance a large number of metallic blades (10 mm 
by 150 mm by 0.1 mm) covered with black soot. The blade is 
then eroded during a fraction of the avalanche time (typically 
1 minute). Taking a photograph of the image of an homoge- 
neous light, reflected by the light, the erosion intensity can 
be determined. As shown in controlled lab experiments, the 
grayscale is a function of the relative velocity between the flow 
and the blade, (b) Picture of the blade after a H — 45 mm 
thick avalanche. The left limit between black and white zones 
correspond to that determined by the string method. The 
grayscale profile is consistent with a plug flow separated from 
the static part of the dune by a shear band. 
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FIG. 6: Power spectrum density (PSD) measured during ho- 
mogeneous steady avalanches similar to spontaneous ones. 
The power radiated is the integral of the PSD. Dotted line: 
avalanche on a 4 m high dune in Tarfaya (d ~ 183 /im). 
Dashed line: H = 65 mm thick avalanche on a 42 m high 
dune in Sidi-Aghfinir (d ~ 163 /im). Solid line: H = 15 mm 
avalanche on the same 42 m high dune. Dotted-dashed line: 
same avalanche (H — 15 mm), the signal coming from an 
accelerometer at 60 mm below the surface. 
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FIG. 7: (a) Schematic of the measurement set-up used to 
characterise surface wave dispersion relation and in particu- 
lar, the cut-off frequency (first resonant mode). The emitter 
is excited with an amplified sinusoidal source. It is composed 
of a coil glued around a sand-paper cylinder buried in the 
sand and submitted to the action of a fixed permanent loud- 
speaker magnet. This way, the force acting on the sand pa- 
per is imposed and, as inertia is small, transmitted to the 
grains. The transducers are piezoelectric accelerometers of 
diameter 10 mm. This size is chosen to have a hierarchy of 
well-separated length scale: the wavelength (typically 40 cm) 
is much larger than the transducer size, which is much larger 
than the grain size d. (b-c) Same as fig. 1 of this comment, but 
on the slip face of a booming dune. Field experiment: propa- 
gation of a gaussian wave-packet at 400 Hz along the slip face 
of a 42 m high mega-barchan. (b) Signal envelope <^(t) mea- 
sured by correlation with the emitted signal at 10 cm (dotted 
line), 220 cm (solid line) from the source, (b) Comparison 
with the model proposed Bonneau et al. 2007. The multiple 
wave-packets received are direct evidences of the existence of 
multiple surface modes guided by the gravity induced index 
gradient. 




FIG. 8: (a) The Makhnovist drum experiment: response of 
the booming dune to a normal tap containing constituting a 
broadband excitation. After the propagative modes excited 
have left, the resonant mode i.e. which does not propagate 
stays. The tail following the tap contains a well defined fre- 
quency that can be interpreted as the first compression reso- 
nant mode, (b) Auto-correlation function of the signal shown 
in (a). The resonant frequency is around 73 Hz for these par- 
ticular conditions of field measurements. The wet sand layer 
is at a depth around H — 50 cm below the surface, [after 
Bonneau et al. 2006] (c) Amplitude £ of vibration at 5 m 
from the source as a function of the frequency /, for perma- 
nent sinusoidal signals emitted at the surface of a 4m high 
booming dune. The wet sand layer is at a depth H — 50 cm 
below the surface. 
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FIG. 9: http://www.pmmh.espci.fr/fr/morphodynamique/Booming.mov 

Movie of homogeneous and inhomogeneous booming 
avalanches in the Atlantic Sahara. From 0" to 38": homo- 
geneous steady avalanche on a 42 m high barchan (d=163 
microns). The sound level in dB can be read directly. 
The graduated stick gives the height of the free surface. 
From 38" to 56" : homogeneous steady avalanche on a 
4 m high barchan (d=183 microns). From 56" to 70": 
homogeneous steady avalanche on the 42 m high barchan, 
recorded by the slider. The frequency is that of sponta- 
neous events as the avalanche is driven by gravity. From 
70" to 88": inhomogeneous unsteady avalanche produced 
by a pulsed slide. The frequency of the squeaks is not 
that of spontaneous events as the avalanche is driven 
by the pressure exerted by the bottom of the slider. This 
movie is to be compared to that provided by the Caltech team 
(http : //www. its . caltech. edu/ nmvriend/research/mult imedia.html) 
that never show a sliding technique able to induce an homo- 
geneous steady avalanche. Contrarily to their statements, 
they do not show either that booming sustains after the 
avalanche stops (one sees the sliders, not the avalanche). 
From 88" to 105": comparison of the sliding techniques. 
[Citation of Auxiliary Material (Animations SI, S2, S3 and 
S4) provided as scientific complements to their article by 
Vriend et al.] 



